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INTRODUCTION 
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i . 

The  NAVSTAR/GPS  NDS-2  satellite  will  be  spin  stabilized 
at  approximately  IDO  revolutions  per  minute  during  its  early 
flight  phases,  but  must  be  nearly  motionless  relative  to  its 
orbit  and  pointing  toward  the  earth  during  its  operational 

I 

lifetime.  Present  plans  call  for  the  despin  to  be  accom- 
plished in  several  phases  allowing  the  nutation  damper  to 
remove  any  nutational  motion  caused  by  thruster  misalignment 
during  the  despin  periods.  This  report  covers  portions  of 
the  USAFA  Independent  Stability  and  Control  Analysis  performed 
during  June  and  July  1976  on  the  dynamic  behavior  of  the 
satellite  during  this  spindown  prior  to  earth  acquisition, 
i Most  results  presented  are  analyt ic,  wi th  results  of  digital 

simulation  included  to  verify  the  analytic  work.  The  last 
two  sections  discuss  simulated  behavior  not  predicted  by  the 
analytic  results,  and  give  heuristic  explanations  for  this 

4 

behavior.  Some  detailed  derivations  and  simulation  summaries 
»*  are  given  in  appendices. 


r' 


II.  BACKGROUND,  TERMINOLOGY,  AND  ASSUMPTIONS 

The  behavior  of  a non- symmetr ic  body  spinning  about  the 
z axis  can  be  represented  by  Euler's  equations.  Using  prin 
cipal  body  axes  these  equations  are 


I U + f I - I ) w w = M 

x x y z yJ  y z x 

1 uj  + (I  - I ) w = M„ 

y y 1 x z x z y 

I to  + ( I ■ I ) w u)  = M 

2 z y xJ  x y z 


(2-1) 

(2-2) 


(2-3) 


u>x,  0)^  and  are  the  respective  inertial  angular  rates  in  the 
principal  axis  frame;  Ix>  I and  Iz  are  moments  of  inertia  about 
the  principal  axes;  M , M , M are  respective  body  axis  torques. 
The  principal  angle  of  interest  in  this  study,  nutation  angle, 
is  defined  in  terms  of  principal  spin  axis  and  tangential 
angular  momenta  as 


...  . HT(t)  * (I/V 

tan  = HTTty  “ 1 a.  ■ 


(2-4) 


where  the  principal  spin  axis  is  the  z body  axis.  In  all 
derivations,  we  assume  v(t)  is  a small  angle,  i.e. 

HT(t) 

tan  v(t)  = v(t)  - — ppy 

z ^ 1 

The  explicit  use  of  t as  an  argument  of  functions  is 


(2-5) 


dropped  where  possible  in  the  interests  of  brevity.  Vve  also 


assume  that  Ix  and  I are  nearly  equal,  i.e, 

d - r ) = ai«i 


(2-6) 


In  addition,  the  small  angle  assumption  for  v implies 

!x  Xy 

IXtV<IZW-  anC*  1 ywy<<  * zWz  ’ l0r  T~  ~ I ~ S0  wx<<(117 

' z 1 z 

and  a)  <<uj  . 

y 2 

The  previous  two  assumptions,  together  with  (2-5),  allow 
us  to  infer  that  in  the  absence  of  torque  about  the  z axis, 

M . 


0 or  a)  (t)  = constant  = co 

7 J / 


(2-7) 


Looking  again  at  the  unforced  case,  we  see  that  if  I „ is 

the  maximum  inertia,  and  we  are  spinning  in  the  positive  sense 

about  the  z axis,  then 

(2-8) 
(2-9) 


co  >0  -=^>  oo  >0 

x y 


co  > 0 =>  to  <0 

y x 


or,  as  shown  graphically  in  Figure  2-1,  we  expect  the 
following  precession  of  the  tangential  (non-spin)  angular 
momentum. 


Figure  2-1 


Natural  Precession  of  Tangential  Angular  Momentum 


This  simple  and  perhaps  elementary  fact  is  quite  useful  in 
interpreting  and  specifying  consistent  initial  conditions 
for  simulations. 
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III.  NATURAL  BEHAVIOR  OF  NUTATION  ANGLE 


In  terms  of  body  angular  rates  and  using  the  previous 
assumptions,  we  may  derive  the  equation  of  motion  for  trans- 
verse angular  momentum  and  nutation  angle  (see  Appendix  A); 


H 

u - _z 
mt  ht 


H 


wxu,yAI  + FT  Mx  + if  My 


• Wxwy  AT  + 1 

v ~Fl^  AI  ht 


H 


H 

+ 


y 


u2  M 


(3-1) 


(3-2; 


In  general,  these  equations  are  difficult  to  integrate 
analytically,  but  some  special  cases  are  enlightening.  In 
particular,  we  first  examine  the  unforced  natural  response 
of  the  system.  By  letting  all  torques  be  zero,  we  are 
neglecting  the  effect  of  a nutation  damper  as  well  as  jet 
torques.  Equations  2-1  and  2-2  have  been  solved  by  various 
authors 

in  most  general  form 


under  these  assumptions  with  the  result  that 


cu  (t)  = A cos  to  t + B sin  o>  t 
y v 1 y n y n 


co  (t)  = A cos  to  t + B sin  <o  t 
x v J x n x n 


where 


(I. 


I Hi 

xJ  y z 


V 


Vy 


co  = K,  co 
z 1 z 


(3-5) 

(3-4) 

(3-5) 


For  GPS/NDS  with  solar  panels  in,  during  spindown,  K1  = .13. 

We  expect  then  the  basic  motion  of  to  and  ui  to  have  a period 

x y 


about  7.6  times  the  period  of  rotation.  Simulation  runs  con- 
firmed this  type  behavior. 

We  may  obtain  relatively  simple  expressions  for  to  (t) 
and  to  (t),  by  defining  the  initial  time  of  our  simulation 

y 

and  t = 0 of  our  analytic  expressions  to  correspond  to  the 

point  on  Figure  2-1  at  which  (0^(0)  = 0 and  to  (0)  = to  = 

7 hnax 


to  . Letting 
7 o 

H,r  = H (0)  = I m 
o 77 

we  have  (see  Appendix  A) 


where 


(3-6) 


(3-7) 


(5-8) 


(3-9) 


We  note  that  if  1^  = I we  will  have  = 1.  K?  is  approximately 

0.7  during  the  spindown  phase  for  GPS,  in  which  case  we  note 

that  to  oscillates  between  smaller  maximum  values  than  to,,, 
x y 

This  behavior  was  also  verified  by  the  simulation.  If  we  now 
put  these  expressions  in  Equation  3-2,  we  obtain  (see  Appendix 


A) 


l 


where 


/T~  + K^(cos  2 - 1) 


(3  in  J 


and 


v(t)  = vq 


v = v (0) 
o 


(3-11  ) 
(3-12) 


(Note  that  if  AT  = 0,  v(t)  will  remain  constant,  i.e.  spin  is 
about  the  principal  spin  axis.)  Simulations  verified  the 
results  of  Equation  3-10.  We  may  also  note  that  vminimum 

occurs  when  id  t = 90°,  270°,  ...,  at  which  time 


J — l,3,o,... 


(3-14) 


These  computations  are  useful  in  that  they  allow  us  to  isolate 
the  behavior  of  v(t)  due  only  to  the  asymmetry  of  the  body,  and 
to  focus  on  those  effects  due  to  torques  such  as  the  nutation 
damper,  thruster  misalignments,  and  intentional  despin  torques. 


TV.  FORCED  BEHAVIOR  OF  NUTATION  ANGLE 


The  integration  of  Equation  3-2,  including  the  effects  of 
forcing  torques,  is  extremely  difficult,  even  under  the  simpli- 
fying assumptions  already  made.  We  may  make  only  general 
observations.  Let  us  predict  system  behavior  as  AI  ■+■  0,  (we 
approach  a body  of  revolution.)  Under  this  most  stringent 
constraint  we  can  show  that  if  we  are  despinning  about  the  z 
axis  and  have  no  torques  other  than  about  the  z axis,  we  expect 


v(t)  = 


s-  v. 


(4-1) 


This  behavior  is  due  only  to  the  fact  that  H7(t)  is  being 
reduced  while  H^,(t)  is  not  (under  these  assumptions).  This 
model  describes  the  "mean”  behavior  of  v(t)  during  spindown 
surprisingly  well,  but  tends  to  overestimate  the  growth  of 
v (t ) as  time  progresses.  It  is,  however,  useful  as  a rough 
guide.  One  effect  for  which  this  mode  does  not  account  is 
that,  for  1. 1 f 0,  at  the  time  despin  firing  stops,  the  nuta- 
tion angle  (as  presented  in  Equation  3-10)  may  not  be  at 
its  natural  maximum.  As  a result,  we  in  fact  notice  v(t) 
increasing  even  after  despin  firing  stops,  sometimes  up  to 
approximately  1/K-,  times  its  value  at  the  time  of  thrust 
termination . 

Another  effect  quite  noticeable  in  simulations  is  the 
sensitivity  of  final  nutation  angle  to  the  time  of  firing 
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when  error  torques  (such  as  thrust  misalignment)  are  present. 


This  sensitivity  should  be  anticipated  from  Equation  3-2.  he 

note  that  M and  M may  increase  or  decrease  the  nutation 
x y 

angle,  depending  on  the  signs  of  the  torques  and  the  phase  of 
the  precession  (signs  of  and  Hy) . Since  the  signs  and 
magnitudes  of  M and  My  are  by  definition  unknown  random  con- 
stants (see  Appendix  BJ,  it  is  not  possible  to  choose  the  "best" 
time  to  fire  for  despin,  even  if  precession  phase  were  measur- 
able, which  it  is  not  for  GPS/NDS.  Since  the  forced  behavior 
predicted  in  Equation  4-1  is  dependent  on  vQ  at  time  of  firing, 
and  this  in  turn  is  dependent  on  the  unmeasured  precession 
phase  angle  during  unforced  rotation,  we  have  an  additional 
random  element  introduced.  With  additional  instrumentation 
such  as  rate  gyros,  it  would  theoretically  be  possible  to 
estimate  the  thruster  misalignments,  as  they  are  unknown  but 


constant . 


Two  effects  could  be  seen  during  despin.  One  increase  in 
v was  a nearly  linear  function  of  v . Equations  4-1  and  3-10 
summarize  this  effect.  The  other  increase  in  v was  independent 
of  vq  and  was  the  effect  of  misalignment  torques  and  My. 

The  increase  in  v during  a thrust  period  was  nearly  linear  in 
its  reaction  to  these  torques.  The  two  effects  cannot  always 
be  directly  added  because  they  may  not  be  "in  phase." 


V.  NUMERICAL  RESULTS  OF  SPINDOWN  SIMULATIONS 


Extensive  simulation  runs  of  spindown  operations  from  20 
rpm  down  to  1 rpm  were  made  with  various  thrust  misalignments. 
No  major  problems  were  encountered  down  to  6 rpm,  as  the  nuta- 
tion damper  is  reasonably  effective  down  to  this  speed,  and 
any  nutation  angle  built  up  is  dissipated  before  starting 
the  next  step  of  spindown. 

The  most  comprehensive  runs  involved  the  region  below 
6 rpm.  Error  torques  simulated  corresponded  to  misalignment 
angles  of  0,  1,  and  2 degrees.  As  expected,  the  larger  error 
torques  gave  the  more  extreme  behavior.  Simulations  were 
made  starting  with  .5°  nutation  angle  at  6 rpm  and  despinning 
to  2 rpm  - then  starting  the  spindown  to  1 rpm  at  several 
different  times  in  the  precessional  cycle,  starting  from  the 
natural  motion  existing  at  the  end  of  spindown  to  2 rpm. 
Another  set  of  runs  with  the  same  range  of  thrust  misalign- 
ment errors  was  made  spinning  down  directly  from  6 to  1 rpm. 
This  will  be  referred  to  as  the  one-step  approach.  Although 
not  a complete  sampling  of  all  firing  times  and  error  mis- 
alignments, Figure  5-1  summarizes  the  range  of  results 
obtained : 


1 1 


; 


Final  Range 
2-Step 
Method 


Final  Range 
1-Step 
Method 


Spin  Rate  [ r pm j 
Figure  5-  1 

Result?  cf  Low-speed  Spindown  Simulation 


The  exact  numbers  for  final  nutation  angle  at  various 
speeds  might  be  different  with  a statistically  significant 
sample,  but  the  general  shape  and  relative  positions  of  the 
final  nutation  angle  ranges  should  remain  the  same. 


1 2 
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VI.  RECOMMENDATIONS  AND  CONCLUSIONS 


Simulation  results,  based  on  thruster  torques  due  to 
misalignment  of  2°  or  less,  indicate  that  most  nutational 
behavior  will  be  of  the  type  discussed  under  natural  behavior 
and  behavior  due  to  intentional  despin  torque.  At  any  spin 
rate,  but  noticeably  at  low  speeds,  the  effect  of  misalign- 
ment torques  may  be  either  stabilizing  or  destabilizing  and 
is  unpredictable  in  both  magnitude  and  sign  unless  further 
instrumentation  is  added  to  the  present  system.  This  should 
be  considered  for  future  missions. 

The  following  conclusions  are  drawn  from  simulations 
summarized  in  Section  V and  specifically  address  the  ques- 
tion of  whether  it  is  advantageous  to  stop  the  despin  process 
at  2 rpm  and  then  proceed  to  1 rpm.  Firstly,  in  a probabi- 
listic (mean  value)  sense,  there  is  no  advantage  to  stopping 
at  2 rpm,  assuming  that  the  nutation  damper  has  negligible 
effect  at  this  speed.  The  time  constant  of  the  damper  effect 
is  on  the  order  of  several  hours  theoretically.  Secondly, 
the  worst  and  best  nutation  angle  resulting  from  a plan 
stopping  at  2 rpm  would  be  larger  and  smaller  respectively 
than  if  despin  was  continued  directly  from  6 to  1 rpm.  If 
the  management  criterion  is  to  minimize  the  maximum  nutation 
angle  rather  than  the  expected  angle,  it  would  appear  better 
to  despin  directly  from  6 to  1 rpm.  Under  the  assumption 
that  we  have  no  nutation  at  6 rpm,  we  should  still  expect  2° 

IS 


i 


..4u  * 


to  3°  nutation  at  1 rpm,  no  matter  how  we  despin,  if  we  have 
approximately  thruster  misalignment.  This  tange  should 
be  recognized  as  an  engineering  judgement  since  the  cost  of 
simulating  sufficient  combinations  of  firing  times  and  mis- 
alignment angles  to  make  a valid  statistical  statement  is 
prohibitive. 
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APPENDIX  A 

DERIVATION  OF  DYNAMIC  EQUATIONS 
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A 

L.' 


Given: 


HT  = 


«L_  (-[ 

dt  T 


_('x“xT  * Vy“yf 

■k[ 


V“x“x  * Iy2“y“yJ 


Using  Euler's  equations  2-1  through  2-3 


— H = — 
dt  HT  H 


T 


m u uj  I A I + u M I + w M I 
xyzz  xxx  y y y j 


To  derive  the  nutation  angle  equation,  we  assume  that  w, 
is  constant  except  for  the  change  caused  hy  Mz . 


e v 


d I HT 


U7Hj  * ^z^T 


d t d t \ H . 


dv 

1 

dt 

“ HzHT 

dv 

i 

e 

i* 

dt 

ht 

!“xVdziI)  * HxMx  * HyMy 


“A 

H 2 
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H, 


v, 
V + A 


H M 
x x 


' Ht  \2 

if  K 

\ z / 


(A-l) 

(A-2) 


(A-  3 ) 


(A-4) 


(A- 5) 


(A-  6) 


Starting  from  the  result  that  for  u>,  constant  and  no  torques,  *x 
and  a)v  behave  as  a second  order  system  with  natural  frequency 

(D  = K,W  (Reference  3).  We  have 
n 12 

u (t)  = Ay  cos  unt  + By  sin  «nt 


(A-7) 


' ~Jt  : '*  ■*  «ABH  . 


WyU)  = -Ay%  sin  V + By“n  C0S  V 

(D  (t)  = A cos  CD  r + B sin  id  t 
x x n x n 

<Dx(t)  = -Axu)n  Sin  unt  + Bx<Dn  cos  u^t 
We  define  the  time  t = 0 such  that 


^y  CO)  = <Dy 


CD  (0)  = 0 

X 

Which  implies  from  the  Euler  equations 

«y(0)  = 0 

Substituting  A- 11  and  A- 13  in  A- . and  A- 8 gives 

B = 0 
y 


Substituting  A-12  and  A-14  in  A-9  and  A-10  gives 
A = 0 


Using  the  definition  of 


B 


x 


I I 
JLX 


(I  -I  )(I  -Ivi 


- CD. 


I (I  -I 

_XLJ— 

icryr 


i 


Therefo re , 


W (t)  = <0  COS  (0  t 

y y n 

1 ’ o 


<o  ft)  = -oo  K-,  sm  oo  t 
x v J y 2 n 

' o 


where  is  defined  by  Equation  A- 19. 


(A- 20) 
(A- 21) 


With  this  convenient  form  for  and  u , it  is  relatively 
straightforward  to  integrate  the  unforced  special  cases  of 
Equations  A- 3 and  A-6.  Since  they  are  essentially  the  same 
equations,  we  will  integrate  only  Equation  A-6. 

dv  _ - A I 


tt  to  2 K-  cos  to  t sin  to  t 

dt  H.p  yQ  2 n n 


-AI  to  2 K, 


y0 


2 Ha 


sin  2oo  t 

n 


(A- 22 ) 


(A- 23) 


t 

#'* 


2to  sin  2oo  t 
n n 


The  bracketed  term  involving  inertias  may  be  reduced  to 
the  simple  form 


K _ A I / IzK2  \ _ 1 
3 ~ 2 ( T 2 I 2 [I. 

y Ki 


AI 


Fx / Vy, 


(A- 24) 


(A-25) 


(A- 26 ) 


(A- 27) 
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so  that 


Intcg  rating 


where 


g i-ving 


vdv  = — j—  K„  sin  2w  (2w  t)dt 
2 o n v n v 

we  obtain 


2 V ' 

v _ o 


K-.  cos  2w  t + C 


n o 


:o  = -f-  C1  ' K3> 


v(t)  = v /T  + K (cos  2 u.1  t - 1) 
o v 3 n 
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APPENDIX  B 


DERIVATION  OF  ERROR  TORQUES 


This  section  develops  the  methods  used  to  determine 
error  torques.  We  assume  the  use  of  two  opposed  thrusters 
(nominally  12  and  14)  which  provide  torque  about  the  nega- 
tive z axis  to  effect  spindown.  They  are  positioned  with 
respect  to  the  center  of  mass  as  shown  in  Figure  B 1. 


Looking  in  detail  at  thruster  12,  we  can  define  misalignment 
angles  in  the  x,y  and  y,z  planes  as  shown  in  Figures  B-2  and 


F 

- - 

I 

z * 

Figure  B-2 
Definition  of  e 

> z 


p 


x 


Figure  B-5 
Definition  of  e 


Jl.J.  * • 


I 


■ 


are  defined  so  as  to  provide  positive 


The  angles  e and  e 
y L A / 

torque?  about  the  x and  y axis  respectively.  We  shall  assume 
that  the  misalignment  angles  allow  use  of  the  "small  angle 
approximation."  We  know  that 

= r x F (B-l) 

where  M and  r are  both  with  respect  to  the  body  center  of 
mass.  We  may  compute  for  thruster  12 


M = 


M 


M 


M 


( B-  2 ) 


ri2  = Th 


xy 

-1 

yz 


where  Th  is  the  thruster  level . Then 


(B-3) 


M = Th 


re  + r 
y yz  z 


r e 
z xy 


r c 
x yz 


- r - r e r 
x y xy 


(B-4) 


We  may  do  the  same  for  thruster  14.  To  keep  the  compu- 
tations more  simple,  we  assume  the  two  thrusters  are  located 
symmetrically  with  respect  to  the  center  of  mass,  but  this  is 
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not  necessarily  true. 


Tf  the  thruster  block  were  rotated,  it 


We  note  that  the  error  torques  are  uncertain  in  both  direction 
and  magnitude,  varying  'with  the  thruster  misalignment  angles, 
which  may  be  independent  of  each  other,  and  possibly  with 
uncertainty  in  thruster  position.  The  misalignments  could  be 
estimated  if  accurate  measurements  of  body  angular  rates  were 
available.  When  Section  V speaks  of  "2°  misalignment,"  it 
means  that  both  thrusters  have  a positive  value  of  2°  in  both 
errors . 
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